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= (cosa/cosP)p + (sina/cos(3)r + tanpcosji y + [siny + tan(3 sin \i cosy] % 
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ROBUST, NONLINEAR, HIGH ANGLE-OF-ATTACK CONTROL DESIGN 
FOR A SUPERMANEUVERABLE VEHICLE 
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ABSTRACT 

High angle-of-attack flight control laws are developed for a 
supermaneuverable fighter aircraft. The methods of dynamic inversion and 
structured singular value synthesis are combined into an approach which addresses 
both the nonlinearity and robustness problems of flight at extreme operating 
conditions. The primary purpose of the dynamic inversion control elements is to 
linearize the vehicle response across the flight envelope. Structured singular value 
synthesis is used to design a dynamic controller which provides robust tracking to 
pilot commands. The resulting control system achieves desired flying qualities and 
guarantees a large margin of robustness to uncertainties for high angle-of-attack 
flight conditions. The results of linear simulation and structured singular value 
stability analysis are presented to demonstrate satisfaction of the design criteria. 
High fidelity nonlinear simulation results show that the combined dynamic 
inversion/structured singular value synthesis control law achieves a high level of 
performance in a realistic environment. 


t Stability and Control Engineer 

* Aerospace Engineer 

* Aerospace Engineer 
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INTRODUCTION 


Supermaneuverablity is defined as the ability to maneuver an aircraft up to 
and beyond the stall angle-of-attack. Some tactical payoffs of high angle-of-attack 
maneuvering include superior survivability, confusion of adversary pilots, and the 
ability to increase first-shot opportunities 1. Additional control power in the form of 
forebody vortex flow control and thrust vectoring can allow fighter aircraft to 
operate in the post-stall flight regime. Therefore, advanced control law design 
techniques must be found for robust high angle-of-attack stability augmentation 
and maneuvering. The resulting controllers must provide for the integration of both 
conventional and unconventional control effectors. 

Modem robust multivariable design methods provide an efficient means of 
developing linear controllers for aircraft. Since the flight control problem is 
inherently multivariable, and the linear aircraft model has associated uncertainties, 
robust multivariable methods are a good choice for flight control design when 
nonlinearities are not too severe. In a Flight Dynamics Directorate contracted 
effort 2 - 3 , a robust H M controller within an inner/outer loop framework was designed 
for a supermaneuverable aircraft at a single flight condition, and robust 
performance was demonstrated for a Herbst-like maneuver. A robust controller for 
this same vehicle was designed by Sparks 4 for a single flight condition using p- 
synthesis in a model-following framework to simultaneously incorporate flying 
qualities specifications and account for structured uncertainty. A recent Wright 
Laboratory technical report 5 describes the design of a p-synthesis controller for a 
supermaneuverable vehicle that is integrated into an inner/outer loop control 
structure to provide full-envelope robust stability and performance for angles-of- 
attack up to 25 degrees. 

Traditionally, flight control law development for low to moderate angle-of- 
attack flight regimes has been accomplished using linear design methods on 
linearized models of the aircraft. However, the advantages of supermaneuverability 
dictate that future air combat will venture into high angle-of-attack, nonlinear 
flight regions. Purely linear controllers are not able to effectively control 
supermaneuverable aircraft for more than very limited flight envelopes. This 
limitation has motivated a number of researchers to explore nonlinear techniques 
such as dynamic inversion. Bugajski and Enns 6 have used nonlinear dynamic 
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inversion to control the HARV aircraft across a wide, high angle-of-attack flight 
envelope. Huang 7 has used a dynamic inversion approach to develop high angle-of- 
attack control laws for the X-29. 

High angle-of-attack maneuvering is still a relatively new area in flight 
controls. Venturing into the regions of post-stall flight can and should elicit serious 
questions about safety issues such as control effector saturation and departure 
susceptibility. Different methods have been successfully demonstrated that assist 
in preventing the destabilizing effects of control saturations. In an approach used 
by Bugajski and Enns®, loop bandwidths are reduced so that a scaled projection of 
the desired control vector is achieved and the control surfaces lie on the boundary of 
an achievable subspace. A method of allocating control effectors such that the 
maximum possible moment is generated within a constrained set of achievable 
values has been suggested by Durham®. Another approach introduces thrust 
vectoring controls when saturations occur in aerodynamic surfaces 5 . 

The main contribution of the work presented in this paper is the integration 
of some the most promising approaches described above into a detailed design 
approach for achieving robust high angle-of-attack flight control designs. The most 
notable advancement is the integration of dynamic inversion and structured 
singular value synthesis. Linearization of the vehicle dynamics is accomplished 
through a nonlinear dynamic inversion scheme. A robust compensator is designed 
around the linearized plant using p-synthesis in a model-following framework. The 
p-synthesis design satisfies flying qualities requirements and robustness goals 
throughout the design etivelope. A control allocation scheme is used which uses the 
pseudo-inverse of the control distribution matrix to allocate controls based on body 
axis rotational acceleration commands. A method known as daisy-chaining is used 
to generate thrust vectoring commands when aerodynamic control effector 
saturation occurs. Adverse control power saturation effects are minimized by 
scaling lateral commands based on an achievable control vector. Control effector 
prioritization is implemented through a daisy-chain technique that limits lateral 
control power demands that compete with longitudinal power requirements. 

In the following sections, a description of a modified F-I8 aircraft model is 
given followed by the definition of design requirements. A brief theoretical 
background is presented on nonlinear dynamic inversion and p-synthesis. The 
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controller architecture and design is described followed by the control allocation 
scheme and departure resistance logic. Finally, linear robustness analysis results 
and the results of a high fidelity nonlinear simulation of a supermaneuver are 
presented. 

AIRCRAFT MODEL 

The aircraft model described in this paper is based upon a modified version of 
the F-18 aircraft. The vehicle is a twin engine fighter aircraft with a moderately 
swept wing, twin canted vertical tails, and a large leading edge root extension. The 
aircraft model is augmented with two dimensional thrust vectoring nozzles that 
provide pitch and yaw moments when deflected symmetrically and a roll moment 
when deflected asymmetrically. The aerodynamic control inputs to the aircraft 
dynamics are the elevators, the ailerons, the rudders, and the leading and trailing 
edge flaps. The aerodynamic surfaces are useful at normal flight conditions, where 
there is adequate aerodynamic control surface effectiveness. The thrust vectoring 
inputs are useful at high angle-of-attack, low dynamic pressure operating 
conditions, where the traditional aerodynamic control effectiveness is inadequate. 
The pilot inputs include a control stick and rudder pedals. 

A nonlinear simulation model of this aircraft exists as modular FORTRAN 
code. The model consists of separate modules describing the atmosphere, nonlinear 
equations of motion, aerodynamics, engines, thrust vectoring nozzles, variable 
geometry inlets, sensors, and actuators which include rate and position limits. The 
high-fidelity model was developed as part of a previous effort which gives more 
detail than that presented here®. There are five pairs of aerodynamic surfaces: 
three pairs for active control and two pairs scheduled for optimum performance. 
The ailerons, rudders, and elevators are used for stability augmentation and flight 
path manipulation. The leading and trailing edge flaps are scheduled to maximize 
airframe performance across the flight envelope. The aerodynamic data are 
contained in tabular format and linear interpolation is used for traditional force and 
moment aerodynamic coefficient build-up. Thrust vectoring-induced aerodynamic 
effects are added to static and dynamic baseline aerodynamic coefficients to obtain 
total aerodynamic coefficients. 
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The dynamics of this vehicle can be described by the following set of first 
order nonlinear differential equations 6 . 10 . The first three equations describe the 
rotational dynamics of the aircraft in the body axis. 


q * ^[™ ne ro + m thrust + P^Z ’ + W** 2 - P 2 > 1 


^aero ^ ^thrust + PQ^xz + qr(Iy * Iz) 


r * i 


lx *Ixz 

-1 

p 






. ‘Ixz Iz . 

- 


( 1 ) 


(2) 


q, p, and r are the body axis pitch, roll, and yaw rates respectively. I x , Iy, Iz» an d 
I xz are the moments of inertia. The m, l, and n terms are the aerodynamic arid 
thrust moment contributions to the rotational equations of motion. 


The next three equations describe the evolution of aircraft motion with 
respect to its velocity vector. 


. „ cosp . . cosy • 

a - q * (pcosa + rsina)tanp - — y- sm\i — x 


$ = psina - rcosa - sinp y + cosp cosy x 


P = 


C0S<X p + ^^r + tanP cospy + [ siny+ tan P sinp cosy] x 


( 3 ) 

( 4 ) 

( 5 ) 


cosp cosp 

a, P, and 1 1 are the angle-of-attack, sideslip angle, and roll angle about the velocity 
vector. 

The last three equations describe the orientation of the velocity vector with 
respect to inertial space. 


- cosp T cosp sinp 
mV 


cosy 


Y = Tiiv L ' mV Y ' V g + mV Fth ' usl Y 


X = 


1 [ sinp L + cosP cosp Y] + —y Fth. «st x 


mV cosy 


V = - “ - g siny + sinpY + ^F lhrustv 


( 6 ) 

(7) 

( 8 ) 
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y is the flight path angle, x is the ground track angle, and V is the velocity. The 
Fthrust terms are the linear contributions of vehicle thrust to the aircraft equations 
of motion , resolved into the respective vectors. L, D, and Y are lift, drag, and side 
force, respectively. The parameter m is vehicle mass, and g is gravitational 
acceleration. 


DESIGN REQUIREMENTS 

Flying qualities are the primary measures of performance for a manual flight 
control system. For conventional flight, specifications for flying qualities can be 
found in MIL-STD-1797A 11 . While these requirements are not valid for high angle- 
of-attack flight, enough guidance is given to provide a basis for extrapolation. A 
rigorous study of new flying qualities measures is beyond the scope of this paper. 
The following high-a requirements are defined here only as baselines for this design 
study. 


Requirements for the short period mode include constraints on the frequency 
and damping of a low order fit of the transfer function between pilot inputs and 
aircraft pitch response. Appropriate forms for this low order transfer function and 
methods for deriving the low order fit are described in MIL-STD-1797A. Short 
period frequency, co sp , should be a function of equivalent airspeed, V^. An 
appropriate guideline for short period frequency is: 

to sp (rad/s) = 1.0 x (9) 

Therefore, the desired pitch response speed to pilot commands should increase with 
equivalent airspeed. At flight conditions above 30 degrees angle-of-attack, short 
period damping should be at the high end of the military standard’s leveT 1 and 2 
requirements. For these conditions, the acceptable range for £ sp is between 0.7 
and 2.0 . 

The primary roll subsidence mode flying quality parameter is roll mode time 
constant, Tr. The roll mode time constant is found from a first order fit of the 
transfer function between pilot input and roll rate response. Past experience with 
fighter aircraft has shown that desired values for Tr are a function of angle-of- 
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attack, MIL-STD- 1797 A level 2 requirements are used as a baseline range of 
acceptable roll mode time constants. The target values used in this study are: T R = 
0.30 at a = 0 deg, T R = 0.75 at a = 30 deg, and T R = 1.40 at a = 60 deg. 

The desired directional response to pilot inputs can be derived from 
requirements on the Dutch roll mode. Dutch roll frequency, coj), and damping, 
can be derived from a second order fit of the transfer function between pilot input 
and sideslip response. Because of the danger of departure susceptibility at high 
angles-of-attack, the Dutch roll damping is required to be greater than 0.8. The 
Dutch roll frequency requirement is taken directly form MIL-STD- 1797 A, coj) > 1.0 

rad/s. 


DYNAMIC INVERSION 


- The purpose of dynamic inversion is to develop a feedback control law that 
linearizes the plant response to commands. In general the nonlinear aircraft 
dynamics can take the form 

i = «x,u), y = Cx dO) 


where x is an n-dimensional state vector, u is an m -dimensional input vector, C is a 
pxn matrix, and y is a p-dimensional vector of output variables. A transformation 
is necessary to put the equations in a form from which the inverse dynamics can be 
constructed. Each controlled output, yj, is differentiated until an input term from u 
appears 12 . Only m outputs can be controlled independently by the m available 
inputs, therefore p must equal m. As shown by Lane and Stengel 12 , the output 
equations may now be written in the form, 


y (dl 


ryiWin 

y 2 ld 2 l 


= h(x) + G(x)u 


L y p Mpl J 


( 11 ) 


where yjldjl represents the djth derivative of the output yj. The inverse dynamics 
control law can be written as 

u = G(x) H V - h(x) ) (12) 
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h(x) represents the nonlinear output dynamics and G(x) represents the nonlinear 
control distribution. The parameter v represents the desired linear dynamics of the 
closed loop system. With the inverse dynamics control law implemented, the closed 
loop system now has the form, 

ytd] ss v 

p 

If the system is observable and Zd^n, then all of the closed loop poles may be 
J i*l 

placed If Z df < n, then closed loop stability cannot be proven. In this case the 
i- 1 ... , 

unobserved dynamics or the internal dynamics of dynamic inversion must be 

checked at local operating points to insure stability 

STRUCTURED SINGULAR VALUE SYNTHESIS 

The structured singular value (p) framework provides a unifying measure which 
can be used to simultaneously address stability and performance robustness 
specifications 1 4 • 1 5 . If p is less than unity for a properly scaled system, then the 
specifications are met. It is desirable to be able to address these multiple objectives 
directly within a design method, p-synthesis provides for the direct incorporation of 
robust stability and performance goals into a design by combining H„, design with 

structured singular value analysis 16 - 17 . The p-synthesis problem is described by the 
attempt to find a controller that minimizes an upper bound on the structured 
singular value, 

min inf sup (j(DM(K)D' 1 ). 

K DeD a 

M(K) is the weighted closed loop transfer function shown in Fig. 1. 

One approach to this problem is the DK-iteration; it calls for alternately 

minimizing sup aOMOOD 1 ) for either K or D while holding the other constant. 
First the controller synthesis problem is solved using H*, design on the nominal 
design model, G w . p-analysis is then performed on the closed loop transfer function 
M(K), producing values of the D scaling matrices at each frequency. The resulting 
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frequency response data are fit with an invertible, stable, minimum phase transfer 
function which becomes part of the nominal synthesis structure. With D fixed, the 
controller synthesis problem is again solved by performing an Hoo design on the 

augmented system. The DK-iterations are continued until a satisfactory controller 
is found or a minimum is reached. Fig. 2 shows a flow diagram for the DK-iteration. 
The resulting controller order is the order Of the design plant and weighting 
matrices, in addition to the order of the D-scale transfer function fits. With each 
iteration, the D-scale frequency response data from the previous iteration is 
combined with the current values, and then the transfer function fit is performed on 
the combined data. This approach avoids a built-in increase in controller order that 
would result if at each iteration new D-scale fit transfer functions were augmented 
into the synthesis model from the previous step. It is important to note that the 
DK-iteration is not guaranteed to converge to a global minimum, but practical 
experience has shown that the method works well for a broad class of problems 17 . 

CONTROL LAW DEVELOPMENT 

Dynamic inversion and structured singular value synthesis are combined to 
achieve robust manual control for high angle-of-attack flight. The controller 
structure is shown in Fig. 3. Pilot pitch stick inputs command a pitch rate, q re f, roll 

stick inputs command a stability axis roll rate,|i re f, and pedal inputs command a 
sideslip, P re f. The following sections describe each element of the control design. 

Fast Inversion 

The aerodynamic and thrust induced moments in eqs. (1) and (2) determine 
the classical linear stability and control characteristics of the aircraft. These terms 
may be expanded into derivative form. 

m aero = l C m« Va + 2 Cmc l q + C ' n &E 5e ^ 2 pVSc 


/aero = t C,p V P + \ C lp p + \ C [r T + C,^ 8 D T 


rtaero = t C„p V P + 2 ^ n P P + 2 ^ nr r + ^' 1 8 dt ^ DT 




+ ^'*6a + ^Sr ^ ^ 

m thrust 

= t C m6pTV Spry! 

(18) 

^thrust : 

- [ C, 5rtv 5rtv + C, 6ytv 5ytv ] 

(19) 

^thrust 

= [C„ 8rtv 5rtv + C n5yTV 5YTV] 

(20) 


where S is the wing area, c is the mean aerodynamic chord, and b is the wing span. 
8 e is the symmetric elevator position, 5 dt is the asymmetric elevator position, is 
the aileron position, 5r is the rudder position, 5prv is the symmetric pitch thrust 
vectoring nozzle position, 8 rtv is the asymmetric pitch (roll) thrust vectoring nozzle 
position, and SyTV is the yaw thrust vectoring nozzle position. The derivatives in 
eqs. (15-20) can be represented in dimensional form where: 


pVSc 2 

“ 4Iy Cfrt 4 ’ 

pV 2 Sc 
M a 8 = 2I y 

L* - pVSt2 C; 
^P.r ~ 4I X w P> r ' 

r • pV2sb r, 
L P,S = 2I X L/ P- 5 

M . pVSb 2 

N P,r = 4I Z L,l P.> ’ 

* pV 2 Sb „ 
%5= 2I Z Cn P-5 


L P,p,r,5 = (Lp.p.r.g + N P,P,r,5 > (J z I x . I xz 2) 


Np.p.r.S = ( N p,p,r,8 +'^ L P,P,r,8 } (1^1^) 


Z 
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( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 
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(27) 
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The rotational equations of motion can now be written as a combination of linear 
and nonlinear contributions. 
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(28) 


At this point, it is assumed that a generalized control scheme has been implemented 
in the form of a control selector, described later in this paper. Eq. (28) can be 
rewritten in terms of the generalized controls: roll, pitch, and yaw acceleration 
commands. 
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(29) 


The first step in control law development is the implementation of a dynamic 
inversion loop that replaces the existing rotational aircraft dynamics with some set 
of desired dynamics. This step is called the fast inversion. 
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Eq. (30) is the application of the dynamic inversion step in eq. (12) to the outputs p, 
q, and r. The stability derivatives and inertial properties in eq. (30) are found 
through linear interpolation of values stored in a tabular database. The desired 
dynamics contain a set of linear stability derivatives that provide satisfactory modal 
frequency and damping characteristics. 
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Structured Singular Value Synthesis 

The fast inversion control law provides equalization of the dominant 
dynamics across the flight envelope. This equalization effectively eliminates the 
requirement for gain scheduling. A robust controller may now be designed around 
the linearizing fast inversion loop to provide command tracking performance. 
Robust tracking of body axis rotational rate commands is achieved with a 
p-synthesis controller. The structured singular value of two uncertainty and one 
performance block is minimized using a DK-iteration. As shown in Fig. 4, three 
frequency dependent weights are chosen to balance performance and robustness 
considerations. 

The design plant is defined by G* = C*(sl - A*) _1 B*. A* is the system 
dynamics matrix at some flight condition which is considered to be central to the 
flight envelope in terms of modal frequencies and damping. B* is the normalized 
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control effectiveness matrix in eq. (29). C* produces the outputs a, P, p, q, and r. 
The fast inversion block represents the body axis rate inversions shown in eq. (30). 

W act defines the uncertainty at the actuator input, W s defines the 
uncertainty at the sensor output, and W p weights the error between the 
complementary sensitivity function of the closed loop system and an ideal model of 
the system response. The actuator and sensor uncertainty models are taken from 
Haiges, et al. 18 . Wp is chosen such that the closed loop system follows the ideal 
model closely at frequencies below 10 rad/s. The ideal model represents the desired 
transfer function between body axis rate commands and roll, pitch, and yaw rate 
responses. For this problem it is defined as 


Ideal 

Model 


~ JL 


" 3 

0 

0 

Pc 


s+3 

9c 

= 

0 

3 

s+3 

0 

r 


0 

0 

3 

- r c - 


s+3 — 


( 32 ) 


A successful p-synthesis design will achieve this first order tracking response to 
body axis rotational rate commands. The diagonal structure of the ideal model will 
also force the response to be decoupled in roll, pitch, and yaw. 

Slow Inversion 

Because only three generalized controls are available, the first step in 
dynamic inversion ignores the dynamics associated with angle-of-attack and 
sideslip. These internal dynamics can be accounted for in a second application of 
dynamic inversion to these slower state dynamics. A simple unitary transformation 
can be made to translate stability axis rate commands into the body axis rate 
commands that are available to the p-synthesis controller. 
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Notice that q c is equal to a c . Eqs. (3) and (4) can be rewritten with these stability 
axis rate commands. 


. . x „ cosn • - cosy 

a = q - (pcosa + r sina ) tanp - 777 ^ 7 - sin M 777 ^X + «c 


cosP 


cosP ‘ 


P = psina - rcosa - sinpy + cospcosyx + Pc 


(34) 

(35) 


If sufficient frequency separation exists between the body axis rate command 
responses and a contribution to the a and P equations, then that contribution can be 

canceled by a slow inversion loop. The inversion control law for a c includes only the 
nonlinear effect of gravitational acceleration due to vehicle orientation. The other 
terms in the a equation are either considered negligible or too fast to be controlled. 


a c = <Jc 


cos[i cosy 
V cosP 1 


(36) 


The inversion control law for the sideslip equation includes the nonlinear term 
representing gravity induced sideslip due to non-zero roll angle. Again all other 
Jterms are considered negligible or too fast to control. 

P c = Kp ( P ref - P ) - Sinn cosy y (37) 


The addition of a sideslip feedback term to this equation provides sideslip command 
tracking and increased turn coordination. The gain Kp is selected to provide a 
second order response that satisfies the frequency and damping requirements for 
the Dutch roll mode. Assume that the closed loop system exactly matches the 
desired first order response in eq. (32). Then it can be assumed that 
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and 
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" s (s + 3) 


(38) 


i . -2_ 

. (s + 3) 
Pc 


P. 

Pc 


With the addition of the sideslip feedback gain Kp, the simplified transfer function 
from P c to P can be written as 


JL 

Pref 


Kp 3 

s (s+3) 3_Kfl 

„ Kr 3 . “ s 2 + 3s + 3K« 


(39) 


A value of 0.5 is selected for Kp, so the sideslip response to commands should be 
second order with a frequency of 1.2 rad/s and a damping of 1.2. 

Command Shaping 

The desired flying qualities for the pitch and roll axis are achieved through 
the use of prefilters. By scheduling these prefilters, the response to pilot inputs can 
be shaped appropriately with flight condition. As described earlier, the desired roll 
response is first order with a time constant that is a function of angle-of-attack. 
With the |i- synthesis compensator implemented, we can assume that the stability 

axis roll rate transfer function is: 


iL . (40) 

. (s + 3) 

Me 


so the response with a first order prefilter with gain is: 



Mref Me 


3Ku 

(s + 3)(s + K^) 


3Ku 

(s 2 + (3+^)s + 3K^) 


(41) 


A schedule for the roll prefilter gain that achieves the desired equivalent system 
response with angle-of-attack is: 

» 3.65 - 0.0433 a and min(F^) = 0.5. (42) 
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The desired pitch response to pilot inputs is second order with a short period 
frequency that is a function of equivalent airspeed. The transfer function 
representing the pitch response to pilot copimands is. 


_0_ 3 (43) 

q c " (s + 3) 


The pitch response with a first order prefilter with gain Kq is: 


w - ^ i 


3K n 


3K n 


(s + 3) (s + K q ) " (s 2 + (3+K q )s + 3K q ) 


(44) 


The prefilter gain, K q , is scheduled to provide the desired level of damping and 
increase in short period frequency with equivalent airspeed. 



(45) 


CONTROL LIMITING AND PRIORITIZATION 

The control selector, sometimes referred to as pseudo-controls, has two functions. 
The first is to normalize control effectiveness by transforming generalized rotational 
rate commands into actuator position commands. The second is to take advantage 
of available control redundancy by allowing for control redistribution without 
changing the linear closed loop performance. The basic idea of the control selector is 
in redefining the control contribution to the state equation (28,29), 

B5 = B*8* (46) 

B and 8 are the actual control effectiveness matrix and control vector. B* and 8* are 
the generalized control effectiveness mat rix and control vector. The actual control 
can now be defined in terms of the generalized control, 


The transformation, T, is the control selector. It is defined simply by 
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T = N(BN) # B* 


( 48 ) 


The operation ( )* is a pseudo-inverse and N is a matrix that may be used to 
combine controls or emphasize/de-emphasize a control channel in the case of 
redundant effectors. Because the B matrix in eq. (48) is a function of flight 
condition and aircraft state, the control selector is a function of parameters such as 
Mach number, altitude, angle-of-attack, and engine power level angle. 

The generalized and actual controls for the supermaneuverable vehicle are 
given by 

5e 
5dt 
5a 

5r < 49 > 

SpTV 

5rtv 
-5ytv- 



Consider the following partitioning of the control effector vector as shown in eq. (28) 


5= [fM' 

where 

^aero = 

5e 

5dt 

5a 

_ 5r _ 

> ^tvec — 

5ptv 

5rtv 

LnvecJ 




_5ytv_ 


resulting in 


B 


aero 


0 L§ dt Lg A Lg R 
M 6E 0 0 0 

0 N$ a N$ r 


B 


tvec 


0 LgRTV L§ytV 
Mgprv 0 ® 

0 NgRTV N§ytv . 


(51) 


353 



With the above partitions, eq. (46) is written as: 


[®aero 


®tvec^ 


*8 



' 1 
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0 ' 

u aero 

= B *8* , where 
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0 
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0 

-^tvec- 


0 

0 

1 


( 52 ) 


A daisy-chain method is used to generate thrust vector commands. Thrust 
vectoring is used only when the aerodynamic surfaces are not able to generate the 
necessary forces and moments required for commanded maneuvers. Therefore, the 
computation of aerodynamic control commands is independent of thrust vectoring 
control commands. The control selector is defined by 


C — T 
°aero = 1 aero ° 


^tvec — ^tvec ^ 


(53) 


and 


^aero = N aero (B aero N aero ) # ^tvec N tvec (B tvec N tvec ) (54) 

where N and N f are used to weight the redundant control effectors. Since the 
aero ivec w 

ailerons contribute more to the roll acceleration and the first priority of the 
horizontal tail should be pitch control, the differential horizontal tail command is 
reduced by weighting the command to be a quarter of the other aerodynamic 
commands. There is no redundancy for the thrust vectoring control effectors, and 
thus, the weighting matrices become 
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Computation of the control selector eq. (53) depends on flight condition. 
Therefore, the elements of B aer0 and B tvec are found using linear interpolation of 

stored table values. 

Nonlinear elements, such as position and rate limits, are required to 
implement the daisy-chain. Fig. 5 shows the structure of the nonlinear control 
selector. A limited aerodynamic surface command (8 aer0 ) is generated from a 
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rotational acceleration command (5*) via the aerodynamic control selector (T aer0 ), 

the aerodynamic surface limits, and command scaling logic. An achievable 

* 

aerodynamic rotational acceleration vector (5 avail ) is computed from the limited 
aerodynamic surface command using the control distribution (B aero ). The difference 
of the commanded and achievable rotational acceleration vectors (eg*) is 
transformed to a thrust vector command (8t vec ) using the thrust vector control 
selector (T tvec ) and command prioritization logic. 

The command scaling logic limits the acceleration command in the event of 
control effector saturation. The lateral/directional generalized control command 
that is generated by the control system can be thought of as a vector. This concept 
is illustrated in Fig. 6. When saturation occurs in one axis, the resulting control 
vector loses both the magnitude and the direction of the desired control. By scaling 
the command vector in both axis, an achievable control vector can be realized that 
preserves the direction of the desired command and holds the limiting controls on 
their limits. A block diagram of the command scaling logic is shown in Fig. 7. The 
figure shows that the scaled vector is the product of the commanded vector and the 
minimum ratio of available and commanded acceleration, 8 sca | e d = 8 x sat. The 
scaling parameter, sat, is always less than or equal to unity. It can be argued that 
when saturations occur, control bandwidth is too high. An interpretation can be 
made that this scaling logic acts to reduce the control bandwidth in the event of 
control power saturation 6 . 

The command prioritization logic limits the amount of commanded 
differential pitch (roll) thrust vectoring. By using models of rate and position 
limiters within a daisy-chain, roll thrust vectoring is commanded only when the 
thrust vectoring nozzles are not saturated due to symmetric pitch thrust vectoring 
commands. For commanded rolls at high angles-of-attack, it can be interpreted that 
rolls commands correspond to performance and pitch commands correspond to 
stability. Therefore, the pitch thrust vectoring command, and thus stability, has top 
priority. 
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ROBUSTNESS ANALYSIS 


A range of flight conditions must be selected for the purpose of linear 
analysis. Table 1 describes the seven conditions that span a broad range of Mach 
numbers, altitudes, and angles-of-attack. Linearized models of the vehicle 
dynamics at these points are used for robustness analysis. 

Table 1 Flight Conditions for Linear Analysis 


Flight 

Condition 

Mach 

Altitude 

(ft) 

Angle-of- Attack 
(deg) 

Equivalent Airspeed 
(ft/s) 

1 

0.2 

10,000 

30 

109.7 

2 

0.2 

10,000 

45 

109.7 

3 

0.2 

10,000 

60 

109.7 

4 

0.2 

30,000 

75 

72.09 

5 

0.4 

30,000 

50 

144.2 

6 

0.6 

30,000 

20 

216.3 

7 

0.6 

30,000 

30 

216.3 


The linear analysis models at each of these test conditions include high order 
actuator models, vehicle dynamics, and control elements shown in Fig. 3. 

The robustness of the closed loop system is tested to simultaneous structured 
and unstructured uncertainties. The structured uncertainties consist of 
perturbations in aerodynamic stability and control derivatives. They are shown in 
Table 2. The structured uncertainties are presented in additive form because 
uncertainty percentages can vary greatly with flight condition, especially when the 
nominal value of a parameter approaches zero. The values shown in parentheses in 
Table 2 are the uncertainty percentages at flight condition 1. These are presented to 
show the relative degree of uncertainty in the different coefficients. Uncertainties 
in thrust vectoring are not included because the control distribution logic dictates 
that those effectors are only used at conditions where linear analysis is no longer 
appropriate. The dimensional form of these uncertainties can be found using eqs. 
(21-27). 
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Table 2 Structured Uncertainty Levels 


stability derivatives 

control derivatives 

AC Z „ = 0.150 (20%) 

ACz5E = 0.0100 (21%) 

AC y p = 0.0150 (12%) 

AC y5 R = 0.0016 (10%) 

AC/p = 0.0030 (10%) 

AC/5DT = 0.0015 (2%) 

AC/p = 0.3000 (66%) 

AC/6DA = 0-0010 (2%) 

AC/ r = 0.1000 (31%) 

AC/sr = 0.0010 (40%) 

AC m « = 0.150 (40%) 

ACmsE = 0.0200 (5%) 

AC mq = 6.0000 (86%) 

AC rt SDT = 0.0015 (6%) 

AC n p = 0.0030 (10%) 

ACnSDA = 0.0015 (6%) 

AC np = 0.1250 (140%) 
AC n , = 0.0750 (24%) 

ACn5R = 0.0020 (5%) 


Unstructured uncertainties include uncertain actuator and sensor dynamics. Fig. 8 
shows the levels of multiplicative uncertainty that must be tolerated for each 
actuator and sensor channel. The quantities were derived as part of the work 
presented by Haiges, et al. 18 . The same level of unstructured uncertainty is 
assumed for all of the aerodynamic control effectors. 

The results of structured singular value analysis indicate that the closed loop 
system is robust to the levels of uncertainty considered. Fig. 9 shows the upper 
bounds for the structured singular values at each of the linear test points. The fact 
that these bounds are less than unity at all frequencies provides a sufficient 
condition for robust stability. The peak in the lateral/directional bounds at 2-3 
rad/sec indicates that Dutch roll mode is the most sensitive to plant uncertainties. 

NONLINEAR RESULTS 

In order to test the nonlinear performance of the flight control system, batch 
simulations are run on a high fidelity six degree-of- freedom simulation of the 
supermaneuverable vehicle. A challenging supermaneuver that tests the 
performance of the control laws and the control distribution logic is a very high 
angle-of-attack velocity vector roll. Fig. 10 shows such a maneuver where the 
aircraft is pitched up to 80 degrees angle-of-attack and then rolled 180 degrees 
about the velocity vector. This supermaneuver creates a rapid 180 degrees change 




in heading angle. The solid arrows represent the aircrafts velocity vector. Fig. 11 
shows the time histories for this maneuver. Actuator positions are given in terms of 
left and right tail (&r L , &r R ), left and right aileron (5 AL , 5 AR ), rudder (5 R ), left and 
right pitch thrust vectoring (SpxvL- 5 PTVr)< an( ^ V aw thrust vectoring (Syrv)- The 
left and right convention is used in place of symmetric and asymmetric so that 
control effector saturations are properly represented. 

The 180 degree change in roll and heading angle is achieved by holding a 30 
deg/sec stability axis roll rate command for six seconds. The performance of the 
dynamic inversion/p-synthesis control system is demonstrated by the smooth, well 
damped stability axis roll rate response and the excellent turn coordination at 80 
degrees angle-of-attack. Less than 2 degrees of sideslip is generated during the 
supermaneuver. Notice that all of the aerodynamic surfaces saturate during this 
maneuver, forcing the control distribution, scaling, and prioritization logic to be 
activated/ Command scaling comes into effect due to rate saturations in yaw thrust 
vectoring at the application and removal of the stability axis roll rate command. 
The pitch thrust vectoring prioritization logic is activated when symmetrical 
horizontal tail saturates, causing a requirement for symmetrical pitch thrust 

vectoring. 

CONCLUSIONS 

High angle-of-attack control laws have been developed for a 
supermaneuverable vehicle with thrust vectoring capability. The methods of 
dynamic inversion and structured singular value synthesis are successfully 
integrated into a design approach which achieves desired performance and 
robustness levels. An advanced generalized controls approach is demonstrated for 
the allocation of redundant aerodynamic and thrust vectoring effectors. Command 
scaling and prioritization are implemented to minimize the destabilizing effects of 
saturations during demanding supermaneuvers. The design goals are achieved 
across a broad range of airspeeds, altitudes, and angles-of-attack. High fidelity 
simulations show that the nonlinear aspects of the control laws perform well in a 
highly dynamic, nonlinear environment. 
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Fig. 2 DK-iteration Process 

Richard J. Adams, James M. Buffington, and Siva S. Banda 
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Fig. 3 Closed Loop System 

Richard J. Adams, James M. Buffington, and Siva S. Banda 
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Fig. 4 Design Model for p-synthesis 
Richard J. Adams, James M. Buffington, and Siva S. Banda 
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Fig. 5 Nonlinear Control Selector 
Richard J. Adams, James M. Buffington, and Siva S. Banda 
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turn coordination 



Fig. 6 Lateral/Directional Command Scaling 
Richard J. Adams, James M. Buffington, and Siva S. Banda 
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Fig. 7 Command Scaling Logic 
Richard J. Adams, James M. Buffington, and Siva S. Banda 
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Fig. 8 Magnitudes of Unstructured Uncertainties 
Richard J. Adams, James M. Buffington, and Siva S. Banda 
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Fig. 9 Upper Bounds on Structured Singular Value 
Richard J. Adams, James M. Buffington, and Siva S. Banda 
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Fig. 10 High Angle-of-Attack Supermaneuver 
Richard J. Adams, James M. Buffington, and Siva S. Banda 
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